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Abstract. This paper is devoted to the study of the gravitational instability of a medium permeated by a 
uniform magnetic field along which a circularly polarized Alfven wave propagates. We concentrate on the case 
of perturbations purely transverse to the ambient field by means of direct numerical simulations of the MHD 
equations and of a linear stability analysis performed on a moderate amplitude asymptotic model. The Alfven 
wave provides an extra stabilizing pressure when the scale of perturbations is sufficiently large or small compared 
with the Jeans length Lj. However, there is a band of scales around Lj for which the Alfven wave is found to have 
a destabilizing effect. In particular, when the medium is stable in absence of waves, the gravitational instability 
can develop when the wave amplitude lies in an appropriate range. This effect appears to be a consequence of 
the coupling between Alfven and magnetosonic waves. The prediction based on a WKB approach that the Alfven 
wave pressure tensor is isotropic and thus opposes gravity in all directions is only recovered for large amplitude 
waves for which the coupling between the different MHD modes is negligible. 
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1. Introduction 

Formation of the largest clouds in the interstellar medium 
(ISM) certainly rely on ram pressure of large-scale 
turbulent flows jVazquez-Semadeni et al.~19 95), whereas 
molecular clouds most likely owe their stability to inter- 
nal turbulent motions. It is thus of importance for the 
understanding of the ISM dynamics, to study the gravi- 
tational instability in media stirred by MHD turbulence. 
Since the first stability analysis by Jeans (Jeans 1902), 
there has been a series of studies of the onset of grav- 
itational instabilities in rotating, magnetized or tur- 
bulent media (e.g. IChandrasckh aT" 19611 IToomre"l 964 
Elmcgre en 1991| ). Neutral turbulence was shown to oppose 
gravity by a turbulent pressure if its characteristic scale 
is small enough compared to the instability wavelength 
IjUhand rasckha r" 1951(1 . In comparison with the mecha- 
nism that generates density fluctuations, turbulent pres- 
sure is however a higher order effect l|Sasao 1973|> . Later 
studies predict a reversal of the Jeans' criterium if the tur- 
bulent spectrum is shallow enough (Bonazzola et al. 1987 
|Vazquez-Semadeni Gazol 1995| . Numerical simulations 
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in two dimensions confirm that highly compressible turbu- 
lence can have both a stabilizing effect for long- wavelength 
perturbations and a destabilizing one at small scales 
IJLeorat et al. 1990|l. 

Consideration of MHD turbulence is crucial in the 
study of star formation. As is well known, magnetic fields 
are ubiquitous in the ISM and they are believed to be 
one of the major agents controlling star formation (e.g. 
IShu et al. 1987|l . If the value of the ambient magnetic 
field (supposed to be relatively smooth over the large 
scales) is large enough (sub-critical regime), overall col- 
lapse is hindered and the interstellar matter gathers un- 
der a quasi-static contraction due to ambipolar drift until 
a super-critical core is formed : this is the low-mass star 
formation process. On the contrary, if the magnetic field 
is weak, matter can easily become gravitationally unsta- 
ble and high-mass stars form in a rapid collapsing process 
l|Mouschovias 1 991). In this picture turbulence is consid- 
ered as an additional pressure against gravity. This pres- 
sure could be provided by turbulent magnetic field fluctu- 
ations. 

In the previous picture the interactions between the 
turbulent motions and the magnetic field are however not 
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consistently taken into account and it is legitimate to won- 
der whether this star formation scenario still holds in a 
fully turbulent case. The key point in this bi-modal star 
formation theory is that magnetic fields always oppose 
gravitational instabilities, at least in the linear regime. 
A uniform magnetic field is known to increase the crit- 
ical wavelength for instability by a factor {(3/1 + 
(where [3 is up to a constant, the ratio of thermal over 
magnetic pressure) if the perturbation is exactly perpen- 
dicular to the ambient field, whereas it has no effect in 
every other direction. Support along the magnetic field 
is provided by Alfven waves (AWs). Recent linear sta- 
bility analyses indeed show that AWs provide an extra 
pressure for perturbations parallel to the background field 
and as a result the Jeans length can be arbitrarily in- 
creased for sufficiently large wave amplitude l|Lou 1 996 
Fukuda k Hana wa" 1999(1 . Numerical simulations in a slab 
geometry have confirmed this result for a spectrum of 
AWs ((Gammie k Ostiker 1996). However, higher dimen- 
sional models including decaying MHD turbulence (e.g. 
IStone et al. 19981 IQstriker et al. 1999jl have shown that 
in magnetically supercritical clouds the time for gravita- 
tional collapse does not depend on the initial level of tur- 
bulence. In high resolution three dimensional MHD simu- 
lations, the presence of short wavelength MHD waves ap- 
pears to be able to delay the collapse ( He itsch et al. 20f)T)l 

Other 



studies by Pudritz (1990) 



McKee k Zweibel (1995) use a WKB 



and 
theory 

(Dewar 1970)| to argue that the AW pressure tensor 



is isotropic and that these waves provide an additional 
support in every direction. The phcnomcnological ar- 
gument is based on the observation that the pressure 

pSvSv reduces 



tensor of the waves Pw = — 

SB 2 



<5B(5B 



to Pw = since for a traveling AW, the equations 

of motion imply that <5v = ± ( 4 ^i/ 2 ■ This conclusion is 
however based on the assumption that the only waves 
propagating in the medium are pure AWs. Even though 
AWs are certainly the only kind of MHD waves that can 
propagate over large distances due to their quasi-non- 
dissipative properties (they are transverse waves), as soon 
as they are perturbed, they generate long wavelength 
magnetosonic waves. The latter arc intrinsically coupled 
to AWs and can drastically change the global response of 
AWs under an external compression. 

In this paper, we address the simple question of the 
linear stability of a self-gravitating medium permeated by 
a uniform magnetic field Bo along which a finite ampli- 
tude circularly polarized AW propagates. In particular, we 
study the case in which perturbations are transverse to the 
uniform magnetic field. We first present a numerical test 
showing that for transverse perturbations, the presence of 
a longitudinal AW does not always provide an additional 
support against gravitational collapse. In order to explain 
this result we have done a linear stability analysis. Due to 
the presence of several kinds of instabilities taking place 
at different scales (e.g. decay instability ((Goldstein 1978 
|Derby 1978| )) and the coupling between them, the three di- 



mensional analysis performed on the complete MHD equa- 
tions appears to be cumbersome and physically unclear. 
In fact this problem, which leads to linear equations with 
periodically varying coefficients, requires a Floquet anal- 
ysis. When the perturbations are parallel to the AW, the 
infinite hierarchy of dispersion relations obtained in this 
way decouples into a set of physically equivalent disper- 
sion relations ( |Jayanti fc Hollweg 1993| ). However, when 
perturbations have an arbitrary direction with respect to 
the propagation direction of the AW, the modes remain 
coupled and the hierarchy of dispersion relations needs to 
be truncated. Even in this case the resulting dispersion re- 



lation is very complicated (Vinas k Goldstein (1991) ob- 
tain a 78 order polynomial). For this reason we perform 
a linear stability analysis on a reduced description of the 
MHD equations based on long wavelength, small ampli- 
tude perturbative expansion <(Gazol et al. 1999|) . 

The plan of the paper is as follows. In Section 2 we 
present the numerical results. The reduced system of equa- 
tions used for the theoretical interpretation of these results 
is described in Section 3. Its derivation is presented in the 
Appendix, and its linear stability analysis in Section 4. 
Section 5 is a conclusion. 

2. Numerical Simulations 

In this section we present direct numerical simulations of 
the self-gravitating 2.5D MHD equations 



dtp + V • (pu) = 
p(d t u + u • Vu) = 



= --Vp 7 + (V x b) x b 

7 

-/3J 2 pV$ 



d t b - V x (u x b) = 
V-b = 
V 2 $ = p-l, 



(1) 



(2) 
(3) 
(4) 
(5) 



where all variables have been adimensionalized using as 
velocity unit, the Alfven speed va — So/(47rp ) 1 ^ 2 , with 
Bq and po being the ambient magnetic field, chosen to 
point along the x-axis, and the average density respec- 
tively. In equation Q 7 denotes the polytropic gas con- 
stant, /3 = c 2 s /v 2 A with c s the sound speed, and J de- 
notes the Jeans number, ratio of the characteristic scale 
L to the Jeans length Lj = (ircl/Gpo) 1 / 2 . The simula- 
tions have been performed using a pseudo-spectral code 
with a resolution of 128 2 grid points, except when other- 
wise specified. 

Previous numerical work in more than one space di- 
mension have mainly addressed the influence of AWs on 
the gravitational collapse in a turbulent context. In or- 
der to clearly identify the AW pressure contribution we 
choose here to study the gravitational instability close to 
threshold in presence of a parallel propagating, right-hand 
circularly polarized AW which is an exact solution even in 
the presence of gravity. We focus on the case of purely 
transverse perturbations and thus choose the longitudinal 
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Fig. 1. | ln/jfe(0, 1)| as a function of time for {3 = 0.5, 
L x /2ir = 0.95, L ± /2n = 1.7 and B = 0.3 

size of the computational domain L x to be slightly smaller 
than Lj so that the system is stable in the direction along 
the uniform magnetic field. The critical transverse size of 
the domain is 



In all subsequent simulations we fix 7 = 5/3, J = 1 and 
just vary the transverse scale L±. The AW of amplitude 
B reads b y + ib z = {uj/k)u y + iu z = B Q e i{ - k *- ut \ p = 1 
and u x — 0, where for the nondispersive case w = k. The 
wave number is usually taken as k = 1 and a noise of 
10 in the first five Fourier modes is superimposed on 
the uniform density. 

We first consider the case (3 = 0.5, for which pertur- 
bations at scales L± > L±j = 2n^/3 are gravitationally 
unstable in the absence of waves. Choosing L±/2ir = 1.74 
and L x /2it — 0.948 a series of six simulations varying 
the wave amplitude has been performed. When Bo = 
the growth rate of the gravitational instability is veri- 
fied to agree within 1% with the theoretical value 0.068. 
Increasing Bo by steps of 0.1, the growth rate is first ob- 
served to increase up to 0.107 for Bq = 0.3. It then de- 
creases until Bo = 0.5 for which the system becomes sta- 
ble. This result confirms that a large enough amplitude 
AW successfully stabilizes the medium against the gravita- 
tional instability. However it is surprising that the growth 
rate increases for small wave amplitudes. 

This suggests to study the behavior of the gravitation- 
ally stable system in presence of an AW. Numerical sim- 
ulations are now performed with a resolution of 64 x 128 
grid points for L x /2ir = 0.95 and L±/2ir = 1.7 so that the 
medium is gravitationally stable. The system remains sta- 
ble in presence of the AW as long as the wave amplitude 
does not exceeds 0.16 but becomes unstable in the range of 
amplitudes [0.18, 0.35]. Figure [T] shows the typical tempo- 
ral behavior of | lnpfc(0, 1)|, where pk(k x , k y ) denotes the 
Fourier coefficient of the density, for an unstable simula- 
tion. In Fig.|3the density field fluctuations p— 1 resulting 




Fig. 2. Density field fluctuations for (3 — 0.5, L x /2tt — 
0.95, L ± /2tt = 1.7 and B = 0.3, at t = 288 
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Fig. 3. I In pfc (0, 1)| as a function of time for (3 — 2.0, 
L x /2tt = 0.7, L ± /2tt = 1.21 and B = 0.3 

from the same simulation are displayed at t — 288. In 
addition to the concentration of mass in the transverse 
direction, a weak longitudinal modulation is also visible. 
Another set of simulations with /? = 2, L x /2-k — 0.7 and 
L±/2tt = 1.21, shows the same behavior with respect to 
the gravitational instability but in this case, the growth of 
the gravitational instability, shown in Fig. is oscillatory. 

Similar results are also obtained when choosing the 
AW wave number on the third Fourier mode. However in 
the case (3 < 1, larger wave numbers or smaller values of 
(3 lead to stronger longitudinal modulation resulting from 
the coexistence with the decay instability. 

In the next two sections we present an analytical in- 
terpretation of the previous results. 

3. Reduced description 

As mentioned in Sect. the analytical study of the trans- 
verse gravitational stability of the system is a very com- 
plex problem. Its solution would require numerical simula- 
tions of the same complexity as the ones done in the previ- 
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ous section. For this reason and in order to gain some phys- 
ical insight, we use a reduced description adapted to study 
the nonlinear dynamics of an AW propagating along a 
strong uniform field coupled with a small amplitude MHD 
flow in planes perpendicular to this field. This reduced de- 
scription has been derived in the absence of self-gravity 
by Gazol et al. (1999) (see also Champeaux et al. (1999)) 
for two different regimes depending on the values of (3 
and its derivation is presented in Appendix for the self- 
gravitating case. A verification of this asymptotic descrip- 
tion has been performed in the context of Hall MHD for 
the development of transverse instabilities on a moderate 
amplitude AW IjLaveder et al. 2002b(l . Here we only con- 
sider the case (3 ^ 1, for which the reduced equations have 
been extended to take into account the coupling with mag- 
netosonic waves. The case where (3 is close to the resonance 
j3 = 1 is more delicate as both the AW and the sonic wave 
must be considered on an equal footing and obey nonlin- 
ear equations. The derivation is based on the existence 
of a small parameter e, squared ratio of the wave ampli- 
tude to the value of the strong uniform magnetic field, 
that identifies with the squared Alfvenic Mach number, 
M\ = Uq/v 2 Ai where uq is a typical velocity fluctuation. 
In the frame of reference of the AW, a stretching of co- 
ordinates allows to study non- linear effects on long times, 
r = e 2 i and large scales, (£, 77, C) = {e{x — t), e 3 / 2 y, e 3 / 2 z), 
for both the AW and the transverse MHD flow. The scal- 
ing and the resulting equations are given in Appendix 
lAl The coupling between the AW and the transverse hy- 
drodynamics is accomplished through the introduction of 
mean transverse velocity and magnetic fields which de- 
scribe the slow transverse dynamics and result from the 
average over the direction of the uniform magnetic field. 

Within the framework of this description, where the 
small-scale magnetosonic waves have been filtered, it is 
possible to study the interplay between the AW and the 
large-scale transverse gravito-magneto-acoustic mode in 
isolation from the other small-scale instabilities. The scal- 
ing is chosen in order to ensure that the transverse scales 
are of the same order as the Jeans length. 

The reduced description presented in Appendix in- 
cludes the Hall effect, which causes the AW to become 
dispersive. This effect is probably unimportant in must 
situations of interest although it can become relevant in 
presence of dust IjRudakov 2001|> . 

4. Linear stability analysis 

A plane wave Boe^ k ^~ UT ^ propagating parallel to the 
uniform magnetic field in a homogeneous medium, is an 
exact solution to equations (|A.36jl - l|A.45|) when its fre- 
quency and wave number are related by u + k 2 /2R l = 0. 
The dispersive effect is kept for generality but it does not 
affect the analysis. The wave is perturbed as 



b = B (l + A) exp i(k£ - lot + ■ 



where the amplitude, A = A(r),£,r), and the phase, 
tfi = tp(rj, C, t) modulations arc real and taken independent 



of the longitudinal variable. Such a perturbation assumes 
that the wave remains circularly polarized, a correct as- 
sumption in the context of a linear stability analysis, but 
that ceases to be valid when strong modulations are con- 
sidered QChampeaux et al. 1997b| ). 

It is easy to see from equation l|A.44|l that if the mean 
transverse magnetic field b does not exist initially, as we 
assume, it is not generated by linear perturbations. The 
other fields are written as p — p^ + p^ x \ where the per- 
turbation separates into oscillating p^ 1 ^ and mean pW 
contributions taken in the form 



P 



(i) = p e i(fee- W r+0) + cx 



and 



p {1) =p(y,(, T ), 



respectively. Note that the amplitude of fluctuating per- 
turbations, p, is a complex number and varies only in the 
transverse directions, p = p R (r),£,T) + ip 1 (77, C r ) while 
pW is a real quantity. After linearizing and projecting on 
the first longitudinal Fourier mode we get the following 
linear system 



(7) 



B (id T <j> + d T A) + ikB (u x + y - 

-~d x (f3p + b x + (3j 2 $) = 

d T u + d± (B 2 A + (36 + ((3 + \)b x + /?j 2 <l) = (8) 

u x -S = (9) 

u x - (35 - pb x - [3j 2 $ = (10) 

(11) 



Bo, 



ikb x + ^(d* ± A + id* ± (j)) = 

Ba 



d T 6+^-(d x (u x -5) + dl(u* x -6*)) =0 (12) 
B 



d T u x - -f (d ± (b x + u x ) + 91(6* + K)) = 



d T b x + — {d±_u* + d*^u) 



Bo 
2 



(d ± (b x +u x ) +d* ± (b* x + u* x ) 



k 2 - 



-— $ + {d m + <9 CC )$ = S + b x 



(On 



(13) 



(14) 

(15) 
(16) 



Note that a truncation in Fourier space is also necessary 
in this case. The order e contribution in the equation for 
the mean longitudinal magnetic field b x , is kept in order to 
allow transverse magnetosonic waves to propagate. From 
equations @, COJ), and (0 we have 

S — u x 5 = 0. 
By separating real and imaginary parts of the remaining 



(g) equations and then assuming that the oscillating, p p 



and non-oscillating, p^, perturbations have the form 
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we obtain an algebraic system of equations from which we 
find 







0.10 




= cj> x 


(18) 


0.05 




l + C e ~ 


(19) 


0.00 


— ^.,=(2/5)"" 










C -0.05 
-0.10 


: "<:., •!:.;> : 




(20) 


-0.15 


7 K 1= 0.45 




-0.20 





iBfiCel 2 



1 - 



Pf 
K\ 

Pf_ 
K 2 



$ = - 



2kn 

hi 
K 2 ' 



where K\ 



K 2 



K 



and 



p- 



P + (ky%K'i) 



l + C t 



1-/3 + 



Pf 



k 2 /e- 

Equations and (ffljl lead to 







4fc 



4k x 

X 2 



(21) 
(22) 
(23) 

(24) 

(25) 

(26) 
(27) 

(28) 



with C* e2 = PC el (l - j 2 /(k 2 /e + K\)) + 1. Rescaling the 
uniform magnetic field and the frequency of the pertur- 
bations by a factor e -1 / 2 , and the Alfven wave number 
k by a factor e 1 ' 2 , i.e. returning to the original variables, 
the coefficients C e , C e ±, and C £ 2 become independent of 
e. We denote the unsealed coefficients by C, Ci, and C2, 
respectively. The dispersion relation, also independent of 
e, reads 



n 4 

-Vl 2 K\ 
+ [P + 1 



cfgi 

16fc 2 

x 2 



Pf- 

K 2 



16/c 2 



Sn 2 



2(/3+l 




(29) 



b„ 



Fig. 4. Effect of varying the perturbation wave number 
near critical values K± c \ = (2/5) 1 / 2 and K± C 2 = 0.5 for 
/3 = 0.5, j = 1 and k = 1.526 

It can be rewritten as 



a 2 



K\ P + 1 



R2 



ftt 2 - 



0j Z 
K 7 



Q 2 



(30) 



From equation l|29|) it is possible to check that for C\ > 
0, f2 2 is always real, so that stability is governed by the 
sign of the last term of this equation. When the system is 
stable in absence of the AW, i.e. K± > K±j = j(P/l + 
/3) 1 / 2 , it can be seen that as K± < K±_ c \, where 



K 1 



P3 2 



+ 3/4' 



(31) 



and when Bq exceeds a critical amplitude Bq c , the value 
for which the last term of equation l|29|l vanishes, the sys- 
tem becomes unstable. 

The first term in the right hand side of equation l|3Ufl 
contains the contribution of the uniform magnetic field, 
whereas the second one is the contribution arising from 
the AW. When the system is unstable in absence of the 
AW, there is another critical wave number K± C 2 for which 
the second term vanishes. It reads 



K 2 - 

A _Lc2 — 



Pf 



P + 3/2 



(32) 



Always assuming C\ > 0, for perturbations with wave 
numbers K± smaller than (respectively larger than) K± C 2, 
the AW has a stabilizing (respectively destabilizing) effect. 
Thus there exists a band of perturbation wave numbers 
[K±c2> K±ci] around K±j, for which the presence of the 
AW has a destabilizing effect. In Fig. 0] we display the root 
n 2 corresponding to the gravito-magneto-acoustic branch, 
as a function of the wave amplitude Bq for (5 = 0.5, 
j = 1, k = 1.0526 and five different values of K± around 
Kj_d = (2/5) 1 / 2 and Kj_ c2 = 0.5 (note that for (3 < 1, 
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Fig. 5. Effect of varying the perturbation wave number 
for 0=2, j = 1, = 0.826 and k = 1.4286 

C\ > 0). It is observed that for K±_ > Kx& (dotted 
line) the medium remains stable, whereas for K± < K_\_& 
(dashed line) the medium becomes unstable when the AW 
amplitude is large enough. For K± < K± C 2 (double-dotted 
dashed line), the medium remains unstable but the pres- 
ence of the AW has a stabilizing effect. The destabilizing 
effect for K± c2 < K± < K± c i is even more pronounced 
for larger wave numbers k although it quickly saturates as 
k increases. 

It is now interesting to discuss whether the previous 
results can be interpreted in terms of an AW pressure. In 
the case of quasi-uniform perturbations, i. e. K j_ -C k and 
K± <C the dispersion relation (|30|1 approximates to 



Vl 2 w K\ [ (3 + 1 



K 2 



-BlCiKi. 



(33) 



Assuming a polytropic dependence of the magnetic pres- 



sure on the density in the form |6| /2 



oc pi 



linearizing 



equations and taking perturbations in the form 

l(T7|) it is easy to see that for (3 — > 0, j w = 3/2. This 
is the value obtained by McKee & Zweibel (1995) in the 
same regime. The factor d {d = (1 - (3)' 1 as e -> 0) 
arises from the coupling of the AW with long-wavelength 
magnetosonic waves excited by the perturbations to the 
AW. 

When £1 is real and f2 <C 1, i.e. close to the stability 
transition for C± > 0, the dispersion relation approximates 
to 



n 2 



K 2 ± [(3 + l- 



K 2 



0- 



3 
4 



Pf 



(34) 



The second term on the right hand side of previous equa- 
tion is negative when Kj_ < Kxd and overcomes the first 
term on the right hand side when Bq > Bq c . In this case 
the effect of the AW cannot be modelled by a pressure. 

When C\ < 0, fl 2 can be either real or complex, de- 
pending on the values of Bq, k, K± and j. As an example 



we choose the parameters used in the second simulation 
discussed in section for which C\ = —3.73. In Fig. |SJ 
fi 2 is plotted as a function of Bq, when it is real. It can 
be seen that for Bq = 0.3, f2 2 is complex, leading to the 
oscillatory instability observed in Fig. |2| Note that in ad- 
dition to the gravitational instability, when (3 > 1, AWs 
are unstable with respect to the absolute filamentation in- 
stability ( |Champeaux et al. 1997a| ILaveder et al. 2002afl . 
which leads to the amplification of the transverse magnetic 
field in filaments aligned with the uniform field. 

5. Conclusion 

In this paper we have shown that a circularly polarized 
AW propagating along a uniform magnetic field modi- 
fies the transverse gravitational instability in two different 
ways. When (3 is not too large (Ci > 0), the AW provides 
an extra pressure if the perturbation wave number K± lies 
outside a band [K± C 2, Kxd]- On the other hand, when K±_ 
lies in this band a moderate amplitude AW leads to an in- 
crease of the gravitational instability growth rate or even 
to the destabilization of a medium stable in absence of 
waves. 

The asymptotic model used in this paper does not al- 
low to study the effect of large amplitude waves but it suc- 
cessfully explains direct MHD simulations results for mod- 
erate amplitude waves. In addition these simulation show 
that at large amplitude an AW is able to stabilize a gravi- 
tationally unstable medium as predicted by a theory based 
on a WKB approach (McKee and Zweibel 1995). The re- 
sults obtained in the present paper concerning the desta- 
bilizing effect of moderate amplitude waves arise from the 
coupling between the AW and the transversely propagat- 
ing magneto-acoustic mode. This effect, which becomes 
less relevant for large amplitude AWs, cannot be recov- 
ered by the WKB theory, which treats the different MHD 
modes separately. 

The present paper only addresses the linear regime. 
The question arises whether the instability triggered by 
the AW saturates early or leads to a runaway situation. 
Numerical simulations of the nonlinear regime, which re- 
quire the presence of dissipation, are planned. It could 
also be of interest to perform numerical simulations in the 
regime of weak Alfven wave turbulence, a strongly turbu- 
lent situation being more difficult to analyze. 
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Appendix A: Reduced equations 

The derivation of the reduced equations is here performed 
for a self-gravitating fluid in the Hall-MHD approximation 
where equation Jjyl is replaced by 



a,b-Vx(uxb) = -|vx ( ^(Vxb)xb 

Ri \p 



(A.l) 



Gazol & Passot: Alfven waves effect on the transverse gravitationaf instability 



The last term of this equation, is the Hall term associ- 
ated with the generalized Ohm's law and Rt denotes the 
non dimensional ion-gyromagnetic frequency. When the 
Alfven wavelength is close to the ion inertial length, this 
term becomes relevant and makes the AW dispersive. 

In order to perform the asymptotic expansion we define 
the stretched variables £ = e(x —t), r\ = e 3 ^ 2 y, ( = e 3 / 2 z, 
and t = e 2 t and expand 



p- 


= 1 + epi 4 


e 2 p2 + +e 3 p 3 ■ ■ 




(A.2) 




= eu Xl + e 


2 u X2 + e 3 u X3 + ■ 




(A.3) 


b x 


= 1 + eb Xl 


+ e 2 b X2 + e 3 b X3 


H 


(A.4) 


Uy 


= e 1/2 (u Vl 


+ euy 2 + e 2 u V3 - 


+ •••) 


(A.5) 


u z 


= ^ 2 (u Zl 


+ eu Z2 + e 2 u Z3 - 


f •••) 


(A.6) 


l) 'J 


= e X '\b m 


+ eb V2 + e 2 by 3 + 


...) 


(A.7) 


h 


= t 1/2 (b Zl 


+ eb Z2 + ( 2 b Z3 + 


...) 


(A.8) 






e$i + e 2 $ 2 + • • 


•)• 


(A.9) 



We also define J cx je 3 ^ 2 , ensuring that the typical trans- 
verse scales are comparable to the Jeans length. When 
expanding the Poisson equation in powers of e, it can be 
rewritten as 



d m + d ( 



(A.10) 



where p = p\ + ep 2 and $ = <!>o + £$1- 

We then expand equations (|H).(f2j). HA.l|) and (gj in 

powers of e. At order e 3 / 2 , we have 



c%ui + d(b x = 0, 



(A.ll) 



where the transverse fields are given by b = b y + ib z and 
u = u y + iu z . In order to include a coupling between the 
Alfven waves and the transverse hydrodinamic motions, 
we include a mean contribution, corresponding to an av- 
erage over the £ variable and denoted by an over-line, in 
the transverse components of the velocity and magnetic 
fields 



bi = &i(£,?7,C,t) + &i(»7,C,t), 
ui = ui(£,i7,C)T") +u 1 (7?,£,t). 



(A.12) 
(A.13) 



Equation i|A.lll) implies that the fluctuating parts (de- 
noted by tildes) satisfy 



ui = -bx. 

At order e 2 , we obtain 

-d{pi + d^u x i + d v {u y i + u y i) 

+d((u z i + Ug\) = 
-dgbxi + d v (u y i + v,yi) + 3q{u z1 + u g \) 
/ \bi + h\ 2 ^ 







d. 



-u xl + f3( Pl + + 



d € b x i + d v (b v i +b yl )+ d c (b zl + b zl ) 



= 
0. 



(A.14) 

(A.15) 
(A.16) 

(A.17) 
(A.18) 



Using equation i|A.14fl and separating the mean values and 
fluctuations in equation (|A.16() . we get 

d^b xl + 3 v b yl + d c b zl = (A. 19) 

d v u y i + d c u zl = 0. (A.20) 

The transverse velocity field is thus incompressible to lead- 
ing order. Defining d± = d n + id^, Eqs. I|A.15I) and (|A.18|I 
then rewrite 

d i (~pi+u xl + b xl )=0 (A.21) 

dlh + d ± bl = 0. (A.22) 

At order e 5 / 2 , separating in the equations for bi the 
mean and fluctuating parts, we obtain 



d T h - ^rdx{uxb*i - u*bi) = 
and 

d T bi - <%(6 2 + u 2 ) + (u x i + b x i)d^bi 
+df: (bi(u xl + b xl )) - ud(b x i 



(A.23) 



--^d±(uib\ - u\bi) + hd^u x i 

^dxihK - b$i) + = 0. 

Similarly, the equation for u\ leads to 

d T u x + d± (l3(pi+j 2 3>) + b x i + (^) 



(A.24) 



(A.25) 



and 



d T ui + d(_ ( (pi - u xl - b xl )bi ) +(p\- u x i - b x \)d^bi 

\bi? - {\bi?) 



-d^{u 2 + 6a) - hudl + u*d±;)b - ^(bdl + b*d^)u 



+8 ± \P{pi+] z <5>) + b xl + 



-^(dl'bi - dj>l) - ^(d* ± h - d x b$) = (A.26) 

where (.) denotes the mean value with respect to the 
£ variable. The solvability condition for eqs. 1A.24I) and 
(IA.26H leads to 

drh + U~ + "xi +b x i)h 
+(- y + u x i + b x i)d^bi + ~d$(biu x i) 



~\dxUpX+j^) + bxl + ^ 



\bx\ 



+^bi(d* ± bx - dJl) - ^(blh-hbt) 



2Ri 



%6r = 0, 



(A.27) 
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which generalizes the usual DNLS equation by the pres- the dominant contributions, except for the equation for 
ence of coupling to both longitudinal and transverse mean b x where order e contribution is kept in order to allow 
fields. 



When pushing to the next order the asymptotic expan- 
sion for the longitudinal fields and combining the equa- 
tions arising at successive orders e 2 and e 3 for these fields 
and e 3//2 and e 5 / 2 for the transverse mean fields, we ob- 
tain the self-gravitating version of equations derived by 



magneto-sonic waves to develop as a consequence of cou- 
pling between parallel propagating AW and transverse 
perturbations, we obtain 



Champeaux et al. (1999) which reads: 



d T b + d ( (b(u x + b x - |) + -b(u x + b x ) 

+ (u + b).V&+^-%6 = 

d T S + -^(uz - 8) + d^(pu x ) + ^ (d* ± (b(u x 



-6 = 



d^b x + d v b y + <9 c fr z = 



bb* + b*b 



(A.28) 



d T S- 

d T Ux 



(b(u x — 6)) + uS + bu 3 
~c - (b{b x + u x )) - 



+ c.c. 



(A.36) 
= (A.37) 

(A.38) 
= Q[A.39) 



bb, 



+ c.c. 



+V ■ (uS + hu x 
1 



9rU x df 

€ 



2R 



d^d^b - d ± b* 
\b + b\ 2 \ 



= 



(A.40) 



(A.29) 



d T b x + -(d x u*+d* ± u) 



^ + ( l 3( 1 -l)-l)^-b x (b x + u x -p) 



+a 5 

+V • (uu x - hb x ) = 
d T b x + - (-9 s 6 a + V • u) 



(d* ± (b(u x + b x ))) 
[d* ± (bu x ) - d* ± {ub x ) 

d T u + d ± [/3{p 1 +j 2 ^)+b x 




(A.30) 



i (d* ± b(u x 



+~(udiu + u*d ± u)-^(bdlb-bd_ L b*) = 

d T b + ^d ± (ub* - u*b) = 
d v b y - 



d c b z = 0, 



(A.41) 

(A.42) 

(A.43) 
(A.44) 



b x ) + c.c. ) + V • (ub x 



hit*, 



-—d e (d* ± b-d ± b*) = 



where oscillating and mean contributions have been sepa- 
rated. Previous equations can be used to rewrite the solv- 
( A - 31 )ability condition as 



d T u (%(u + b) 



v[p(p 1+ j 2 <b) + b x + C-±- 




1: 



1 



d T b + d i [-(b + b)P+ (u x + ~b x - -5)b - -d ± P 



+Q-VQ = (Vxb)xb 
d T b- -«%(u + b)- Vx(uXb) = 

dtb x + ^(d±(b*+b*)+d* ± (b + b) 



+ (u + b).V6+— %6 



2 2 




1 



(A.45) 



P 



1 



(2b x + \b + b\ 2 -(\b + b\ 2 )) 



(A.32) 
(A.33) 
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